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Abstract 

Motivated by cosmological applications like electroweak baryogenesis, we develop a field the- 
oretic approach to the computation of particle currents on a space-time dependent and CP- 
violating Higgs background. We consider the Standard Model model with two Higgs doublets 
and CP violation in the scalar sector, and compute both fermionic and Higgs currents by means 
of an expansion in the background fields. We discuss the gauge dependence of the results and 
the renormalization of the current operators, showing that in the limit of local equilibrium, no 
extra renormalization conditions are needed in order to specify the system completely. 
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1. Introduction 

In a recent letter, we have presented our results for the computation of the expectation 
value of particle currents in a two Higgs doublet model, in the case in which the Higgs field 
background is space time dependent and CP- violating [l]]. The original motivation for such a 
computation comes from the so called spontaneous baryogenesis scenario @ at the electroweak 
phase transition. The basic ingredients of this mechanism are CP violation in the Higgs sector, 
which requires at least two Higgs doublets, and thick bubble walls, so that at least some of 
the more rapid processes (top Yukawa and gauge flavour diagonal interactions) can be assumed 
to be in local equilibrium in the background of the electroweak bubble walls and a adiabatic 
approximation can be applied. Of course, also in this scenario one must require that the baryon 
number violating processes go quickly out of equilibrium in the broken symmetry phase, which 
seems by now very hard to achieve, at least in the Standard Model J3J and in the minimal 
super symmetric extension of it || . 

Both in the case of explicit and spontaneous || CP- violation, a space-time dependent 
relative phase 5(x) between the two Higgs fields can be present inside the wall of the expanding 
bubble. Due to the fact that the background is not space-time translational invariant, the phase 
5(x) can give rise locally to nonvanishing expectation values for the particle currents present 
in the theory. If baryon number violation is then in equilibrium, these nonvanishing particle 
currents can rearrange to generate a baryon asymmetry. 

In order to reliably compute the final baryon asymmetry, it is then necessary, as a first 
step, to determine the values of the currents induced by the background. In the approach of 
refs. [|, U a rotation of the fermionic fields is done to make the Yukawa couplings real. As a 
consequence, a derivative coupling of the form 

C int ~ 8^8 J», (1) 

where J M is the current corresponding to the rotation, is induced from the kinetic terms. In the 
original paper only the time derivative 5 was taken into account. 5 acts as an effective chemi- 
cal potential, called "charge potential", and particle densities acquire nonvanishing expectation 
values rii ~ qi 5 T 2 , where is the charge of the i-th particle under the given rotation. The 
thermodynamical evolution of the system in presence of baryon number violation then leads to 
a nonvanishing baryon asymmetry. The gradient part of 5(x) in eq. (1) was subsequently taken 
into account in ref. |§ where the role of particle diffusion was also discussed. Nevertheless, 
as we already pointed our in using the interaction term in eq. (1) as a starting point to 
compute the perturbations to the thermal averages presents some drawbacks. 
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Since the phase S is communicated from the Higgs to the fermion sector through the Yukawa 
interactions, any perturbation in the fermion densities rii should vanish in the limit of zero 
Yukawa couplings hi. Also, they should vanish in the limit of zero vacuum expectation value 
for the Higgs fields H^x) = v i(x) exp[i6i(x)] because no spontaneous CP violation is present 
in the Higgs sector in this limit. Naively, one could then expect a suppression factor of order 
(h 2 v 2 (x)/T 2 ), where hi is the relevant Yukawa coupling, for the perturbations in the fermionic 
particle number with respect to the original result. Since one is interested in regions of the 
bubble wall where sphalerons are still active, i.e for values of Vi(x)/T typically smaller than one, 
then the above mentioned suppression factor might be crucial for the scenario of baryogenesis 
outlined above |J. 

The reason why these suppressions do not appear in the original treatment, is that consid- 
ering eq. (1) as the only effect of the background is equivalent to perturbing around the Higgs 
field configuration S(x) = 0, Vi(x) ^ 0, which is not a solution of the field equations. On the 
other hand, taking the field equations into account, one immediately sees that the expected 
suppression factors are recovered also in the original treatment, so that the result comes out to 
be rotation independent, as it should, also in this context. Indeed, by partially integrating eq. 
(1), we obtain a 'perturbation term' which is given by the hypercharge violating terms in the 
lagrangian, coming from the fermionic mass terms and from the Higgs potential. Also, from the 
field equations one can see that d^5(x) vanishes as Vi(x) 2 for vanishing v,i(x). In other words, 
it is the perturbation (1) itself, if properly considered, to be vanishing in the limit of vanishing 
top Yukawa and Higgs couplings, or in the limit of vanishing values for the background Higgs 
fields. 

In the limit in which some interaction rates (like for example those for the top Yukawa pro- 
cesses) are so high as to be in thermal equilibrium inside the bubble walls, the local equilibrium 
values for the corresponding particle numbers are of course independent on the reaction rates. 

On the other hand, since their values are fixed by that of the perturbation term (1), they 
are still suppressed by the above mentioned factor of (hfv (x) 2 /T 2 ). 

In this paper we want to describe in more detail our approach to calculate the expectation 
values of a composite operator O(z) in a CP- violating Higgs space-time dependent background. 
For definitness, we will work on the background of a bubble wall, as the one described above, 
however the formalism is quite general and can be easily extended to consider other interesting 
situations, like for example a non trivial background also for the gauge fields. 

The method relies on a functional expansion in power series of where the coefficients 

of the expansion are the n-point 1PI Green's functions with one insertion of the operator O(z) 
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computed in the unbroken phase. 

We shall discuss the dependence of the expectation values from the gauge in which the 
background ^(x) is expressed, showing explicitely that the expectation values of gauge invari- 
ant operators are gauge independent. Moreover, the problem of current renormalization will 
be addressed. In general, the computation of Green's functions with insertion of a composite 
operator requires the introduction of new counterterms besides those necessary for the renor- 
malization of the basic lagrangian. This fact leads to the well known phenomenon of the mixing 
among the different renormalized currents of the theory. We shall discuss the mixing matrix 
for the renormalized currents present in the Standard Model model with two Higgs doublet 
and use the non renormalization properties for the conserved currents in order to reduce the 
number of independent counterterms. 

Finally, we will discuss the limit of local equilibrium, showing that, in this case, the only 
expectation values that can be consistently computed are those for the conserved currents, so 
that no new renormalization condition is needed to specify the system in this limit. 

The paper is organized as follows: in sect. 2 our field theory method is described on 
general grounds and particular attention is given to the classical equations of motions and to 
the question of gauge invariance. In sect. 3 the explicit calculations for the fermionic and 
Higgs currents in the specific model under consideration are given in details. They are relevant 
for the spontaneous baryogenesis mechanism and their consequence have been discussed in ref. 
|]J. Sect. 4 deals with the issue of current renormalization, whereas the role of the conserved 
currents in the thermodynamical limit is described in sect. 5. Finally, sect. 6 presents our 
conclusions. 

2. The expansion. 

In this paragraph we will discuss our approach to the computation of the expectation value 
of an operator O(z) on a non-zero background for the fields of the theory. 

Our starting point is the finite temperature generating functional for the 1PI Green's func- 
tions with insertion of an operator 0{z) (in the following 0{z) will represent a particle current) 

T [^(x), A(x)] = W [Ji{x), A(x)] - W ^xJj{x)^{x), (2) 

3 

where ^(x) are the classical fields of the theory and Ji(x) the corresponding sources, while 
A(x) is the source for the operator 0(x). Note that the Legendre transformation has been 
performed only on the fields and not on the operator. 
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The quantity we are interested in is the expectation value of the operator O(z) on the 
background given by the fields ^(x), which we will specify later, 

5T [$?, A] 



(O(z)) 



*?(*) 



6A(z) 



= 0[^(x)](z). (3) 



A=0 

We can expand the functional O [<&?(#)] (z) in a power series of $^ 

n=0 ti,...,in ^" 

where the coefficients of the expansion are the n-point 1PI Green's functions with one insertion 
of the operator 0(z) computed in the unbroken phase 



n (n) , ,_1 5 n+1 T A] 



2 ^ 1 (^ 1 ).--^ n (x n )5A(^) 



(5) 

$c = A = 



Now we have to specify the background fields, which give the starting point for the 

expansion in eq. (f|). A priori, the relevant background is given by the solution of the full field 
equations, 

5T A = 0] 



= (6) 

cj> c =<J> c 



with appropriate boundary conditions. In practice, however, we have to consider some ap- 
proximation of this complete, and unknown, solution. A possibility could be to consider the 
solutions of the classical equations of motion as the zero order approximation. However in many 
interesting cases the radiative corrections to the effective potential, either at zero or at finite 
temperature, are crucial in order to determine the shape of the potential and consequently of 
the background. A typical example is that of a first order phase transition induced by finite 
temperature corrections, in which the existence of bubble solutions is due to the T^O radiative 
corrections which induce a second minima in the effective potential. Following refs. || |IIJ, we 
will then consider the classical equations of motion in which the tree level effective potential is 
replaced by the radiatively corrected one. 

In the following, we will consider the Standard Model with two Higgs doublets, Hi and H 2 . 
Hi has U(l)y charge —1/2, and couples to the 'down-type' right handed fermions, whereas H 2 
(charge 1/2) couples to the 'up- type' ones. Moreover, since only top and bottom quarks will 
be relevant in the following we will not consider the lighter quarks and leptons. 

Neglecting all the Yukawa couplings but the one for the right handed top, h t , the classical 
lagrangian is given by 

C = + {V^V^Hi + {V^H 2 ) ] V^H 2 (7) 
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where the most general tree level scalar potential is given by 



V{H U H 2 ) 



(8) 



V = m 1 2 \H 1 \ 2 + m 2 2 \H 2 \ 2 -(m 3 2 H 1 H 2 + h.c.) + \ 1 \H 1 \ 4 + \ 2 \H 2 \ 4 

+ A 3 |i^i| 2 |^ 2 | 2 + A 4 |iYiiY 2 | 2 + [A 5 (iYiiY 2 ) 2 + Ael^l 2 ^!^ + Ayl^^iJxiJa + h.c] . (9) 

In eq. (|j) we have also omitted tha SU(2) gauge part, which will not be relevant in the following. 
Note that of the two phases of the Higgs fields, 9\ and 9 2 , defined by Hi(x) = Vi(x) exp[19i(x)}, 
only the gauge invariant combination 5 = 9i+9 2 appears in the scalar potential. The orthogonal 
combination, x — @i — @2, is the gauge phase, corresponding to the would be Goldstone boson 
after spontaneous symmetry breaking. One can now write down the equations of motion coming 
from the lagrangian in (|8[), with the tree level scalar potential V(Hi,H 2 ) replaced by the 
(finite temperature) radiatively corrected one, V(Hi, H 2 ). We assume that the parameters of 
the effective potential (and the temperature) are such that it has only one minimum in the 
charged directions, given by H + = H~ = 0, whereas there are two minima in the neutral Higgs 
directions. Then, a bubble solution exists, with v i t2 (x) changing from a non-zero value vf 2 
inside to zero outside, and H = everywhere. 

At zero order, we will also take the fermion fields and currents to be zero, so that the 
background is given by the solution of the following field equations: 



<9V, 2 + - 7^1,2 (<V =F £> M X) 2 
dvio 4 







d^v = - [(vf - v 2 2 W5 - K + v 2 2 Wx 



(„2 + i%ff>8 - {v\ - v 2 2 )V»x 



dV_ 



where the gauge invariant quantity V^x is given by 

T^ixX = d^x + eA^. 

Assuming also that there are no electric or magnetic fields at zero order, i.e. F f 
chosing the unitary gauge = (see below), we obtain from eq. (pi]) 



so that eq. (|T2| ) now reads 



vld^ = vjd^, {A„ = 0) 



(10) 

(ii) 

(12) 

(13) 
0, and 

(14) 
(15) 
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Inserting the solutions of the above equations into the expansion in eq. ([|) we will now be 
able to compute the various contributions to the expectation values of the currents as a loop 
expansion. In the following paragraphs we will see that the Higgs currents get a contribution 
already at the tree level, whereas, consinstently with our assumptions, the contributions to the 
fermionic currents, and to (d^F^), arise only at one loop. 

Before concluding this section, we want to show explicitly that the expectation value for a 
gauge invariant operator, as obtained by the expansion in eq. is independent on the gauge 
in which the classical background has been computed. In our discussion above, we have looked 
for a solution of the field equations with no electromagnetic fields, i.e. with = 0, so that 
it was possible to chose a gauge in which = everywhere (see eqs. ([T^JT5|) ). Chosing a 
different gauge gives rise to a different background, in which is diferent from zero, and we 
want to show that the result for the gauge invariant operator O(z) is the same as in the original 
background. 

Following the usual procedure for the derivation of the Ward identities it is straightforward 
to obtain the relation 



a oa 



i 6 ~ ( - - V° 2 — n + H°*JV\ = (16) 



where a is the gauge parameter, a M , H® 2 , and H® 2 * are the classical gauge and neutral Higgs 
fields, and we have put the charged Higgses and the fermion classical fields to zero. Differenti- 
ating the above expression with respect to the source A(x) keeping the background fields fixed, 
and then putting A(x) = we obtain the relevant identity containing the functional O of eq. 

9 ™ - it ( H «M - H f JO _ flj *° + = o. (17) 

"da,, 2\ l 5H^ 1 5HT 2 6H$ 2 5 Hi* J y 1 

The Green's functions entering the expansion (|J) have to be computed in the 'unbroken 
phase' i.e. for vanishing background fields §\ = {H® 2 , H® 2 *, a^}. Taking = in eq. (|17|) , 
one gets 

0, (18) 



d — 

11 Sa^ 



$ c =o 



while, differentiating with respect to H®, 

e 50 



i 



2 6H?{y) 



5^(x-y) + d, S '° 



'■6a,(x)5H?(y) 



= (19) 

$ c =o 



and similar identities are obtained differentiating with respect to the other fields. 
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Now, let's indicate with $„• 



,2 5 fl M 



{ H l,2, H \, 

c f Q / f* 

tions discussed above, and with & { = {H 1 2 , H 1 2 



0} the 'bubble' solution of the field equa- 
-l/e<9 M #} the solution obtained from 
the former by a gauge transformation. Then, we consider the functional O on the background 



-3E-C' 



<3>j and expand it according to eq. (^). Since a'^ ^ in this case also Green functions involving 
gauge fields in the external legs will contribute to the expansion whereas, if we consider the 
background only the Green functions with neutral Higgses in the external legs contribute. 
However, by means of the above identities, one can easily show that there are cancellations 
among these new contributions, so that the result is the same as for 1^. Indeed, at the linear 
order in a'^ in the expansion we have 



50 



5a IM (x) 



1 



<i> c =o 



SO 



- / rf4 x<9 M 
e J oaAx) 



9{x) 







(20) 



where the last equality is due to eq. 
50 



$ c =o 



-FrO' 



d x 



d x 



At the linear order in H 1 we have 
SO 



5Hl{x) 



H,(x)(l-t6(x)/2) + 0(e 2 



(21) 



The contribution 0(6(x)) is canceled by the H 1 a'^ term in the expansion by virtue of eq. (|T9|), 

5 2 



I 



d xd y 



- Id xd y d 



Sa,(x)SH°(y) 
5 2 



$ c =o 



5 aix {x)5Hl{y) 



Hl(y)9(y) + Ol9 



$ c =o 



(22) 



/ 



d x 



50 



5H?(x) 



H 1 (x)e(x) + o(e 



$ c =0 



The contribution of n-th order in 9 in eq. (^Tj), is canceled by the terms coming from the 
H 1 ( a|J [i — . . . n) orders in the expansion, due to the appropriate Ward identities obtained 
by differentiating eq. (|T9l) with respect to a^. 

This result generalizes to the other terms of the expansion: all the contributions of the same 
order in the primed Higgs fields and of any order in a'^ in the expansion for sum up to 

the corresponding term of the same order in the non-primed Higgs fields in the expansion for 
£>[$•], so that 



(23) 
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Since in practical applications one has to truncate the expansion to a finite order in the 
fields, it is clear from the above discussion that this can be done in a sensible way only if we 
consider the background $ i5 where the classical gauge field is zero, i.e. in the unitary gauge. 

3. Computation of the currents. 

We will now apply the method illustrated in the previous paragraph to the computation 
of (Jt( z ))<s>°(x)i the expectation value of the current operator Jf (x) (i = H® 2 , H ± , t LR , b^n) 
on the background Qj k = {H®, H®* , H®, H®*}, which, for definitness, we assume to be the 
'bubble' solution to the field equation of motion discussed in the previous paragraph. 

In the following, we will assume that there are no chemical potentials. Then, it is easy 
to realize that a necessary condition for the presence of non-vanishing expectation values for 
the currents is a non trivial space-time dependence of the background. Indeed, if H^ 2 (x) were 
constant then the only two vectors available for the construction of ( Jf (z))$c would be z M and 
(the quadrivector which identifies the motion of the thermal reference). By using translational 
invariance, we are left with (Jf (z))<t>c ~ const x u^, but this is forbidden by CPT. 

Also, by applying a charge congiugation transformation both to the operator and to the 
background, we see that 



i.e. the expectation value for the current would vanish also if the background were real. From 
these two considerations one can conclude that the expectation values we are looking for, will 
depend on space-time derivatives of a complex background, so that the two phases 9\{x) and 
9 2 (x) discussed previously, (or 8(x) and x( x ))> constrained by fll4j) , will play a crucial role. 

Computing the expectation value as a functional expansion on the background, as in eq. 
(£D, the problem is reduced to the calculation of 1PI Green's functions with one Jf (x) insertion 
in the unbroken phase, as in eq. fl5|). 

We will work in the imaginary time formalism, so that our expressions for the Green's 
functions and for the currents are valid in the reference frame of the thermal bath, where the 
bubble wall is moving. 

The zero order of the expansion in eq. (5) vanishes since it is given by the 1PI part of 
(J|*(z))$c = o which is zero since translational invariance holds in the unbroken phase. 

The linear order is also zero, as can be realized observing that in the unbroken phase the 



(J? ( Z ))^(x) = -(Ji(z))*c {x y, 



(24) 



theory has still a global U(l) invariance, whereas 
symmetry. 
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H?(x) H°(y)* 



to 



Figure 1: The tree level contribution to the neutral Higgs current. 
The first non- vanishing contributions then come from the quadratic terms i.e 

1 r ^ 3 r 

± X ^ / -.A -.A 



^(x)^(y). (25) 



Using again the global U(l), one can see that the only possible non- vanishing con- 
tributions to the above expression come from the 1PI part of (Ji(z)Hi(x)H2(y))$?=o, or 
{Ji(z)Hi(x)H®{y) )$c=o (i = 1, 2). These Green's functions will now be computed as usual in 
perturbation theory, and the result inserted into the expression (]2"5j). 

We start with the Higgs currents = i (^HjV^Hi — V^HjH^j . In the case of the neutral 
Higgses a contribution already appears at the tree level, see Fig.l. Inserting this contribution 
into the expression (p5|) we get the classical current 

(r H0 (z))^ = -2vf(z)d^(z). (26) 

i 

Note that, as expected, the individual Higgs currents are zero if the phases are constant. 
Moreover, from the equation of motion (0), we see that the total Higgs hypecharge current 
vanishes at the tree level. 

At one loop there is a contribution to the neutral and charged Higgs currents given by 
graphs like that in Fig. 2. 

Since the computation has to be performed in the unbroken phase, we must use resummed 



propagators for the Higgs fields in order to deal with the IR divergencies |TT| . In the unbroken 
phase, the Higgs spectrum contains two complex electrically neutral fields and two charged ones. 
At the tree level, the squared masses of one of the neutral states and of one of the charged ones 
are negative, since the origin of field space becomes a minimum of the effective potential only 
after the inclusion of the finite temperature radiative corrections. The resummation can be 
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H?(x) 



H^(x) 

.0 




Figure 2: The 1-loop contribution to the neutral and charged Higgs currents. 

achieved by considering the propagators for the eigenstates of the thermal mass matrix, which 
has positive eigenvalues given by 



M 2 2(T) = m\{T) + m 2 (T) T #j(T) - m 2 (T)) 2 + 4 mf(T) ^ 

where the m 2 (T) are the thermally corrected mass parameters of the potential @, m\(T) ~ 
m 2 + 3g 2 T 2 /16, m|(T) ~ m 2 . + h 2 T 2 /4, while m 2 (T) receives only logarithmic corrections in T, 
which were computed in ref. ||. Correspondingly, the neutral complex eigenstates are given 
by 

h = cos a Hf + sin a if^ * 

(28) 

H = -sinatf? + cos a #2* 
where 

2m? (T) , . 

sin 2a = — ^7— 3V /, 9/ ^ x . 29 

Completely analogous formulae hold for the charged eigenstates. 
The loop in Fig. 2 gives 

Mp|«M=-' As SiD2a ^^/If^-^to. (30) 

where the function IT'(g) is the analytical continuation to Minkowski external momentum q of 
the finite temperature integral, 

H " W) = PA wtmW)] m-py + mn " ((M ' (T) " M22(T)) ) ' (31) 

where q A = iq , q { = q u p = 27in p T, p 2 =pl + \p\ 2 , and 

d 3 p 



(2tt) j 
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Note that the only quartic coupling contributing to the Higgs currents is A5, whereas there are 
no contributions from A6 and A7. Moreover, the contribution vanishes if there is no mixing 
m\{T) between the two Higgs gauge eigenstates (see eq. (^)). 

It is easy to check that Green function in ([31]) is ultraviolet finite, due to the minus sign 
between the two terms in the integrand. After some algebra, the integral in p can be casted 
into the form 

1 



H"(g) = \<? 



(q 2 + M 2 (T)) fdp 



(p» + M|(T))((g -pf + Mf(T))((q + P y + M 2 (T)) 

(M 2 (T) «-> Ml (T) 



(32) 

Now, the scale of the external momentum q is set by the bubble wall width, which, in the case 
of interest for us of thick bubble walls, is given by L~ l ~ T/100 — T/10 0. Then, since we 
have q 2 <C Mf 2 (T) <C (2nT) 2 , we approximate the integral by neglecting the q dependence in 
the integrand and by considering only the zero Matsubara mode for pq. Inserting the result 
into the quadratic term of the expansion (|]) we obtain 

Each charged Higgs gets a contribution equal to that of the neutral Higgs belonging to the 
same doublet. As expected, also the Higgs currents vanish in the limit of vanishing Vi(z). 

The first contribution to the fermionic currents arises at one loop and is given by graphs like 
that in Fig. 3. Once inserted into the expansion (f|), they give, in the case of the left handed 
top, 

<J&(*)> (3) = ih\ J d*xd%Im(H 2 (x)H° 2 (y)*) 0"(s,y,z), (34) 
where h t is the top Yukawa coupling and 

v ' y ' ; J (27r) 4 (27r) 4 v ' ; K ' 

is the Green function corresponding to the diagram in Fig. 3, with 

q^,D = ^[(^ww-M (36 ) 

V 7 T (k + m) 2 k 2 {k-l) 2 V ; 

and the zero component of the fermionic loop momentum is k° = (2n + l)nT. 

Contrary to the case of the scalar loop, this fermionic loop integral is now infrared finite, 
since the zero component of the loop momentum never vanishes. As a consequence, we will 
work with the tree level fermion propagators, without of resummation. 



11 




Figure 3: The 1-loop contribution to the left-handed top current. 

Again, in presence of a thick wall background we can perform a high temperature ex- 
pansion. This is achieved by approximating G M (m, /) ~ G^(m, 0) + G M (0,Z) (note that 
G M (m, I) = — G M (7,m)), and neglecting terms of order l 2 /T 2 in the computation of G^(0, 1). 

Now the integral has a logarithmic ultraviolet divergence, so renormalization is needed (see 
also the next section). In the MS scheme we obtain 

<^)> (3) - vl(z)d»0 2 (z) (log (j^j + fj , (37) 

where A f = vr 2 exp(3/2 - 2 lE ) ~ 13.944. 

Eq. ( |3"T| ) shows the expected dependence on h 2 and v 2 (z) which, in comparison to the 
original result given in ref. M, gives a suppression factor 0(/i t v 2 /7r T) 2 . 

A graph similar to that in Fig. 1 for the right handed top quark leads to a contribution to 
(Jt R (z)) given by (Jf R (z))^ = —(Jt L (z))^. For the other fermion species, one finds analogous 
results, in which h t is replaced by the appropriate Yukawa coupling, and Vi{z) (vi(z)) appears 
for the up (down)-type fermions. 

A contribution to (Jt L (z)) proportional to Im^H^H^) = t>it>2sin5, and to hf, appears at 
two loops, given by the graph in Fig. 4. The corresponding integral is 

'« = fo v + ^ + M ?( r)) tf + Mf(r)) G " ( -« - k ' < 38) 



where q is the bosonic loop variable and G M (m, I) has been defined in (|36D . Computing the 
fermionic integral in the approximation described above, we are left with a ultraviolet finite 
bosonic integration which receives important contributions only from the zero Matsubara mode. 
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Figure 4: The 2-loop contribution to the left-handed top current. The scalar internal lines can 
be either neutral or charged Higgs fields. 



Inserting the result into the expansion we obtain, again in the high temperature approximation^, 

Tm§(T) / / f \ ^ 
64vr 3 (MriT) + M 2 {T)f \ *>\A f T 2 )^2 j 



<#(*)> W - ~^3 ^^l ml Z\^ k (-^) + I) * MzMz) tin6{z)] . (39) 



As for the one- loop result, the two loop contribution with neutral Higgses in the internal 
lines for (Jt R (z)) is opposite to the one for (Jf* L (z)). This is no longer true when the graphs with 
charged Higgses in the internal lines are taken into account. In such a case, (Jf (z)) gets a a 
contribution equal to eq. (|39"D , whereas the result for {Jt L {z)) is analogous but proportional to 



hi instead of hf. The charged Higgs loops give also rise to a non- vanishing left handed bottom 
density opposite to eq. (p9|). 

When applying our calculations to the interesting case of spontaneous baryogenesis, one has 
to remember that the typical values of Vi(x)/T are smaller than one, as required in regions of 
the bubble wall where sphalerons are still active. This fact enables us to stop the expansion 
series at the second order, being the remaining terms suppressed by powers of Vi(x)/T with 
respect to it. 



1 Note a factor 6 of difference with respect to the result of ref. [Q 
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4. Current renormalization. 



In the previous section we have seen that in order to calculate the finite temperature average 
of the particle currents, the expansion (||) can be used, whose coefficients are the 1PI Green's 
functions with the insertion of current operators Jf , computed in the unbroken phase. With 
Jf 's we have indicated the current for the i-th particle, expressed in terms of the renormalized 
fields. 

In general, the computation of Green's functions with the insertion of a composite operator 
requires the introduction of new counterterms besides those necessary for the renormalization 
of the basic lagrangian (see ref. |l!| for a thorough discussion). For example, the one loop 
computation of the Green function Q^{x y y : z) = ( J/^(z)i^2 ( x )-^2 '*{y) ), defined in eq. (|35|) , 
gives a divergent part (in the MS renormalization scheme, and in the momentum space) 

G**»(ro, = +^ t 2 ^ ~ K - n • (40) 

where m, I are the momenta associated to the external Higgs fields and e = 4 — n. This 
divergence may be expressed as 

g^(x,y,z) = --J-^i (JUz)H° 2 (x)H° 2 *(y)) (41) 

16 7T Z e 2 

it i.e it is proportional to the tree level Green function involving the H® Higgs field current, 
corresponding to Fig. 1. 

This is an example of a quite general phenomenon: a renormalized current [Jf], can be 
defined as a linear combination of the currents Jf 

[Jt}=E( S ij + 5Z ij) J j ( 42 ) 
j 

in such a way that all the possible Green's functions with [Jf] insertions are finite. In the 
example discussed above, we have SZi 2 = — fi 2 / (16tt 2 e). The 5Zij are in general new renor- 
malization constants which are not contained in the counterterm lagrangian, so that new renor- 
malization conditions are needed. 

An important exception to this is the case in which J M is a conserved current for the 
renormalized lagrangian 

£ — £basic + £>cti (43) 

where Cbasic is expressed in terms of renormalized quantities and C ct is the counterterm La- 
grangian. Let's consider here the simple case of a lagrangian containing only left and right 
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handed top quarks and the neutral Higgses, which will be enough for discussing the renormal- 
ization of the Green's functions relevant for this paper. 
The conserved current under the group U(l) is then 

J» + SJ» = q L (l + 5Z L ) r L + q R (1 + 5Z R ) r R + ft (1 + 5Z 1 ) J£ ? + q 2 (1 + 5Z 2 ) J^, (44) 

with q 2 = q R — ft = —ft = 1/2 and is the part of the current coming from C c , t .. The 5Z a 
(a = L, R, 1, 2) are now the wave-function renormalization constants for the different fields, 
coming from £ ct . In this case, being the current conserved, no new counterterms besides those 
already present in 



(and hence in £ ct ) are needed W2 



i.e. 



J" + sr = [r] = q L + q R [jft + ft r H 



•hi 



Q2 



TV 



(45) 



Recalling this property of the conserved current it is possible to reduce the number of linearly 
independent renormalization constants 5Zij (where now i = L, R, 1, 2) necessary in order to 
define the renormalized currents according to eq. fl4"2"p. Indeed, from Q4"5j), recalling (52) and 
(fill), the number of new linearly independent renormalization constants 5Z^s is reduced from 
16 to 8. 

The most general renormalization matrix is then given by 

/ 



1 + sz LL 
SZ L - 5Z LL 

sz 1L 

V 5Z 1L + 5Z LL - 5Z L 



5Z R - 5Z RR 
1 + SZ RR 
5Z 1R 

SZ 1R - 5Z RR + 5Z R 5Z 1 



sz L1 
-sz L1 

1 + 5Z n 

f 6Z n + SZ L1 



SZ L2 \ 

sz L2 

5Z 12 

1 + 5Z 2 + 5Z L2 + 5Z 12 J 



(46) 

We see then that, in general, eight extra renormalization conditions are required in order to 
define the set of renormalized currents 

The only divergent graph encountered in sect. 3, is the fermionic one loop (see eq. (37)), 
so, for the present application, and to our order of computation, it is enough to compute the 
matrix ([I6|) only to order (g)°hf. To this order, only one of the eight extra renormalization 
constants is non zero, SZ L2 and moreover we obtain, again from ([|5|), that SZ L2 = —5Z 2 (+ 
higher order terms), so that no new renormalization condition is needed. 

The mixing matrix is then given, at this order, by 



( \Jl\ \ 

m 
v m j 



( 



\ 



i 

5Z L 


-5Z L 



sz R 
1 


6Z R 



6Z 2 

1 





( J £ \ 










/ 





(47) 



where 6Z L2 = -5Z 2 = (-{i 2 h 2 t /l6 vr 2 e), 6Z L = 5Z R = (-fi 2 h 2 t /32 tt 2 e). 
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Working at higher order, we would need new renormalization conditions, defining the physi- 
cal currents [Jf ] at a certain scale. Due to the interactions, at a different scale the renormalized 
currents will be given by different mixings of the unrenormalized ones, so that, in general, the 
definition of 'pure' (i.e. non-mixed) currents is scale dependent. As we have already stressed, 
this is not the case for the conserved currents, which do not suffer any mixing and do not 
require any new renormalization condition. 

As we will discuss in the following section, the only relevant currents in the limit of local 
thermodynamical equilibrium are just the conserved ones, so that in this case, the state of 
the system is completely determined without the need of introducing extra renormalization 
conditions. 

5. The limit of local equilibrium. 

Up to this point, we have computed the expectation values of the individual particle numbers 
starting from a partition function of the form 

Z[Ji, Aj] = Tr exp 

where H is the Hamiltonian of the system and we are now considering only the zero components 
(Nj = Jj) of the currents considered in the previous paragraphs. We assume that all the 
chemical potentials are zero. 

If we take generic values for the sources Jj and Aj the partition function fl4~8"|) will not corre- 
spond to the one describing the limit of thermodynamical equilibrium, in which the individual 
particle numbers assume values such that the free energy of the system is minimized. In the 
case under consideration of a space-time dependent Higgs field background, the limit of local 
equilibrium is achieved only if all the interaction rates are much higher that the inverse of the 
typical space-time scale of the variation of the background, i.e. the bubble wall width. This 
condition is fulfilled in reality only by a few interactions (the gauge flavour conserving ones, the 
top Yukawa interactions, and possibly some Higgs-Higgs interactions), however it is anyway 
instructive to consider first this limit before discussing the more realistic situation in which 
some of the interactions are out of equilibrium inside the bubble wall. 

The free energy of the system is obtained by performing the Legendre transformation of the 



13 H + 



AjNj 



(48) 
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functional defined in eq. (Q) with respect to the sources Aj of the density operators, 



$ t c , Nf\ = W [Jj, Aj] -Ei/ d 4 zJ,(x)^ c (x) - Ei / ^A,(i)^ 



(49) 



= (H) - TS. 

The last equality in the above equation is straightforwardly derived recalling that the entropy 
S is given by S = —dW/dT and that W = —TlogZ (we are taking the three dimensional 
volume Q = 1) . 

The classical densities Mf(x) are defined as usual, 



sm,A 3 



6Aj(x) 



(50) 



A,=0 



where the background is, as before (eq. (||)), the solution of the field equations with appropriate 



boundary conditions, i.e. we have J{(x) = in eq. (f48|). 

Now, the local equilibrium situation on the background ^(x) is obtained once all the 
processes have driven the free energy to its minimum possible value. This condition corresponds 
to a system of M equations (where M is the number of processes in equilibrium) of the form 



(0 



5F 



$ c J\f c 



]>>j%(x) = (/ = 1,...,M), 



(51) 



where Uj is the multiplicity of the j-th particle in the Z-th process (i.e. in the process A + 2B <-> 
C + D we have z/4 = — vc = —vd = 1, and z/g = 2). 

The system (]5Tf) is of course strictly analogous to the usual system of equations that one 
has to solve to determine the equilibrium chemical potentials. In this case, it restricts the set 
of the linearly independent sources to {A^} (a = 1, . . . , N), corresponding to the set of the N 
conserved charges N a of the system. This means that the only composite operators that appear 
in the equilibrium partition function 



Z eq \J^ A'J = Tr exp [-/? (h + j A' a N a 
are those corresponding to the conserved charges, 



A'=o 



(52) 



x 



5A' a (x) 



(53) 



A'=0 
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where T eq is obtained from Z eq in the standard way. 

After solving the system (|5l"l) , the sources for the individual particle numbers, Aj, are 
expressed in terms of the new sources A^ as 

A^E^AL, (54) 

a 

where q$ is the a-charge value of the particle i. Then, the equilibrium values for the individual 
particle numbers are given by 

a UL - 1 a\ Jy J A' =0 

a 

Concerning renormalization we see that since now the only Green's functions that have to 
be computed are those involving the conserved charges, no new renormalization conditions for 
the composite operators are necessary when one is dealing with the local equilibrium case. 

However in many interesting applications not all the processes have rates fast enough such 
that a complete local equilibrium situation can be attained. In this case one can consider an 
adiabatic approximation, in which the 'slow' processes are freezed out inside the bubble wall, 
while the fast ones are in equilibrium. As a consequence, on the bubble wall only the equations 
corresponding to such 'fast' processess will contribute to the system ( pip and one can repeat 
the above analysis with new effectively conserved charges now playing the role of the iV a 's. 

There is in this case a, mainly conceptual, problem, since these charges are now not con- 
served by the full theory, but only by the effective one obtained by sending to zero all the 
couplings contributing to the 'slow' processes. Then they do in general mix one another under 
renormalization, as was discussed in the previous section, so that the interpretation in terms 
of physical currents will be scale dependent. 

However, since the mixing is due only to the small couplings which are neglected in the 
effective theory, it will be in most of the practical applications a negligible effect. 



6. Conclusions 

Motivated by the recent interest on the possibility of generating the baryon asymmetry in 
the early Universe during the electroweak phase transition and, more in particular, by the so 
called spontaneous baryogenesis mechanism in which the baryon asymmetry is created inside 
thick CP-violating and expanding bubble walls of the Higgs fields, we have developed a general 
method to calculate the expectation value of a composite operator in the presence of a generic 
nonvanishing background for some of the fields of the theory. 
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The method is based on a functional expansion of the expectation value in powers of the 
background fields $i(x) which are a solution of the equation of motion with properly chosen 
boundary conditions. At any order in the functional expansion, the expectation value is deter- 
mined by a 1P1 Green's function in the unbroken phase, which can be computed in perturbation 
theory. 

Using Ward identities, we have shown that the expectation value for a gauge invariant 
operator is independent on the gauge in which the background has been computed. 

The renormalization of the currents has been also discussed, stressing the crucial role played 
by the conserved currents in the case of local equilibrium. In this limit, all the individual 
particle currents are determined in terms of the conserved ones, whose computation do not 
require the introduction of new renormalization constants besides the ones already contained 
in the renormalized lagrangian. 

In practice, the system is well far away from such an idealized situation, with the possible 
exception of the processes mediated by the top Yukawa and the gauge flavour diagonal couplings. 
In this case, a better description of the evolution of particle numbers can be achieved by 
employing a system of kinetic equations. For each particle number, the values expressed by 
(p5p would then represent the infinite time limit of the corresponding rate equation, 

Mi ~ Tij(J\fj - Afj q ). (56) 

As we have already stressed in the introduction, the dependence of such equilibrium values on 
the coupling constants and on the value of the Higgs background should not be confused with 
the dependence of the rates T^- on these quantities. 

The computation of the currents is of course only the first step towards a reliable determina- 
tion of the final value of the baryon asymmetry in this scenario, the major source of uncertainty 
being represented by the rate of the baryon number violating processes inside and in front of 
the bubble walls. 

In any case, due to its generality, we believe that the method illustrated in this paper can 
be easily applied to similar situations of interest in Cosmology, in which some of the fields 
acquire a space-time dependent classical value as for instance in inflationary models or 



in presence of domain walls, as considered in in ref. |14| in the case of the Next-to-Minimal 
Supersymmetric Standard Model. 

The problem of the computation of the currents on a CP violating background has been 



addressed recently also by Huet and Nelson [15fl for the case of a fermionic axial current. 



Although they adopt a semiclassical approach quite different from the one discussed in this 
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paper, their results are in quantitative agreement with ours. In particular the same dependence 
on h t and v(x) is found. 
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